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Bose- Fermi Kondo model with Ising anisotropy: cluster- Monte Carlo approach 
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Ahst.rar.t. 

D The Bose-Fermi Kondo model (BFKM) captures the physics of the destruction of Kondo screening, which is of extensive current 
t , interest to the understanding of quantum critical heavy fermion metals. There are presently limited theoretical methods to study 
"^ the finite temperature properties of the BFKM. Here we provide some of the consistency checks on the cluster- Monte Carlo method, 
; which we have recently applied to the Ising-anisotropic BFKM. We show that the method correctly captures the scaling properties 
j^ of the Kondo phase, as well as those on approach to the Kondo-destroying quantum critical point. We establish that comparable 
Ch Jesuits are obtained when the Kondo couplings are placed at or away from a Toulouse point. 
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• The sub-Ohmic Bose-Fermi Kondo models are of con- 
siderable interest in the context of quantum critical heavy 
fermions [1] and certain mesoscopic structures [2]. The fi- 
pite temperature scaling properties of the BFKM have been 
studied using an e-expansion [3] and also at certain large-N 
limit [4]. At finite N and e, standard Monte-Carlo meth- 
ods [5] can be used for the Ising-anisotropic BFKM but the 
lowest temperature that has been reached is about 0.01 T^, 
where T^ is the bare Kondo scale. Recently, we have ap- 
plied a cluster-Monte Carlo method to this problem [6], 
which is able to reliably reach temperatures of the order of 
).0~^ — 10~^T]^. The purpose of this paper is two-fold. First, 
we demonstrate the consistency of this method for both the 
Kondo phase and the quantum critical regime. Second, we 
address the effect of the deviation from the Toulouse point 
on the finite-temperature scaling properties; previous stud- 
ies of the Ising-anisotropic BFKM [6] have focused on the 
Toulouse point of the Kondo couplings. 

The cluster-Monte Carlo method builds on the well- 
established understanding [7] that the Kondo problem in 
the scaling regime can be studied by a one-dimensional 
Ising model with long-ranged interactions. Through a 
Coulomb-gas picture of spin flips, the Kondo Hamilto- 
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mapped onto a classical Ising chain Hj = ^^ KnnSfS,^ 
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with algebraically decaying (long-ranged) 
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interaction K^ {n — rj ) — K/ 1 »'i — tj p~^ placed at its lower 
critical dimension (e = 0) [7,5]. Here, S is the impurity 
spin and s is the electron spin density at the impurity site 
of the c-electrons with a featureless dispersion in its kinetic 
term, Hq{c); S^ is an Ising variable at the chain site i. The 
coupling constant Knn depends on Jj_ while K is soley a 
function of Ju and can be expressed entirely through the 
scattering phase shift 5 of the electrons, 4 tan (5 = t^J\\P, 
where p is the conduction electron density of states at the 
Fermi energy. An extension of this equivalence has recently 
been used to address the quantum critical properties of 
Ising-anisotropic BFKMs [8,6], 



^bfkm - 'T^K + 9/ ,S' 
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^ Wp (t)^<t>p , (1) 



where the impurity spin S interacts with fermions c|,g. and 
bosons (j)^ . The spectrum of the bosonic bath is taken to 
be sub-Ohmic (0 < e < 1), '^ [5{u} — ojp) — 6{uj + ujp)] ~ 
^sgn(w) and gives rise to an interaction Kf^{ri — Vj) = 



along the Ising chain on top of the Kf^ ( 



9/\' 

rj) = K/\ri — rj\'^ for the Fermi-only Kondo model. 

Before discussing the determination of the critical cou- 
pling (jc or gc, it is worthwhile discussing the g = case 
in greater detail. The RG flow of the Kondo model is to- 
wards an SU(2) invariant fixed point on trajectories with 
Jm(&) — Jj_{b) = c, where b parametrizes the RG flow and 
c is a positive constant. At the Toulouse point, the Kondo 
model can be mapped onto a resonant level model [10]. De- 
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Fig. 1. The scaling function for the Kondo model away from its 
Toulouse point: Jj_ = 0.75ro, (1 - 25 /ir)'^ = 3/5 and tq = 1/16. The 
inset shows our results right at the Toulouse point, {1 — 25/tt)-^ = 1/2; 
all other couplings, and the temperatures, are the same as in the 
main plot. 

viation from the Toulouse point is needed to restore the 
SU(2) symmetry at the fixed point. The Ising chain Hj at 
e = does not possess a continuous symmetry but under- 
goes a Kosterhtz-Thouless-Hke phase transition; the RG 
flow resembles that of the Kondo model away from the 
Toulouse point. This raises the question of whether the uni- 
versal properties of the BFKM, obtained from simulating 
the classical Ising chain, are sensitive to whether the Kondo 
model is placed at or away from its Toulouse point. Previous 
simulations have utilized the Toulouse point, (1 — 26 /tt)'^ ~ 
1/2, where a comparison with exact results is possible [5,6] . 
Here we present results away from the Toulouse point of 
Hbficm in Equation (1) with g — 0. We choose (1 — 26 /n)^ — 
3/5, but keep J± = 0.75ro at the value considered in [6] for 
easy comparison. This results in a bare Kondo temperature 
T'^ — 0.553ro, with a Trotter parameter tq = (3/L — 1/16. 
The Ising chain was solved via a Cluster-MC scheme [6,9]. 
Fig. 1 displays the results for the dynamical susceptibility 
x(t, T) =< TrS''{T)S^{0) > for various temperatures. AU 
curves collapse on one scaling curve W{x) with W{x) ~ x'^ 
for X <C Tfc when plotted as a function of the single vari- 
able X — T/ sin^TTTT). This is a consequence of the con- 
formal invariance of the Kondo Hamiltonian [Equation (1) 
with g = 0]. This result immediately implies the correct 
frequency and temperature dependences of the dynami- 
cal susceptibility. A similar scaling collapse applies at the 
Toulouse point, as seen in the inset. 

The static susceptibility versus temperature T for the 
Kondo model [g — 0, away from Toulouse point with (1 — 
26 /tt)'^ = 3/5] is shown in Fig. 2. These results in conjunc- 
tion with the results from reference [6] establish that the 
universal properties obtained via our method are insensi- 
tive to whether the Kondo model is placed at or away from 
its Toulouse point. We can now study the approach to the 
Kondo-destroying QCP at g = gc- As g increases, the static 
local susceptibility deviates from the Pauli behavior of the 
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Fig. 2. Static susceptibility for tq = 1/16, J± = 0.75to, 
(1 — 25/tt)'^ = 3/5 versus temperature for the Kondo model g = 0, at 
an intermediate coupling gjTj^ = 0.545 and at the critical coupling 
gc/Tx = 0.773 for e = 0.4. The divergence of x(9c) is cut off due to 
the finite system size L = /3tJ~ ; the cutoff temperature decreases as 
TO decreases [6]. 

Kondo phase. As g approaches gc, an intermediate temper- 
ature range opens up, which shows a power-law divergence. 
The critical exponent is consistent with the perturbative 
RG result, x{T) ^ T^'^^'''\ The lower cutoff temperature 
of this critical regime becomes smaller as g increases to- 
wards gc, consistent with the general expectations of quan- 
tum critical scaling, xi'^j T) shows a scaling collapse similar 
to that shown in Fig. 1, providing evidence for an enhanced 
symmetry at g = gc, the quantum critical point [6]. 

In summary, we have established that the cluster-Monte 
Carlo study of the classical Ising chain with long-range in- 
teraction can be used to obtain universal properties of the 
Ising-anisotropic Bose-Fermi Kondo model, both at and 
away from the Toulouse point. Because of the large tem- 
perature range it is able to reliably access, the approach is 
expected to be useful for studying the quantum critical be- 
havior of not only the BFKM itself, but also the extended 
dynamical mean field theory of the Kondo lattice model. 
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